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Abstract: We have analyzed pair production in the innermost region of a two-temperature 
external soft photon comptonized accretion disk. We have shown that, if the viscosity parameter 
is greater than a critical value a c , the solution to the disk equation is double valued: one, 
advection dominated and the other, radiation dominated. When a< a c , the accretion rate has to 
satisfy (th < ltl< ltl 2 in order to have two steady state solutions. It is shown that these critical 
parameters itli, rtl 2 are functions of r, a, and 9 e , and a c is function of r and 9 e . Depending on the 
combination of the parameters, the advection dominated solution may not be physically 
consistent. 



INTRODUCTION 

The idea that pair production could occur in 
astrophysical plasmas is not an old one. It dates from a 
pioneering paper by Bisnovatyi-Kogan, Zel' dovich 
and Sunyaev 1971 [2]. In their study, these authors 
have considered the production of pairs through 
particle-particle interactions and have discovered that, 
under these conditions, there is a maximum 
temperature above which productions overwhelms 
annihilation, equilibrium is no longer possible and 
proton density approaches infinity. Soon, after 
BKZS's[2] conclusions on pairs on relativistic 
plasmas, a lot of works followed (Stoeger 1977 [15]; 
Pozdnyakov, Sobol and Sunyaev 1977 [13]: Liang 
1979[6]) which failed to confirm the existence of this 
maximum temperature. However, though considering 
different processes for pair production, these works 
have focused on a rather limited range of the 
parameter space. A full understanding of BKZS's[2] 
results was only possible after the work of Lightman 
1982[8], who showed the limitations of BKZS's[2] 
results. Accordingly, these results only hold for very 
small photon densities, or, equivalently, very small 
scattering depths. Much before attaining the maximum 
temperature, photon-photon processes will dominate 
pair production. In that regime, the relation between 
pair production and annihilation is quite different, not 
allowing for a critical temperature as is the case of 
particle-particle interaction. As a matter of fact, 
Lightman [8]has shown the existence of two branches 
of solution for a given (x p , T e ): one, the low n+, with 
particle-particle dominating pair production, and the 
other, the high n+, with photon-photon interactions 
dominance in pair production. In the upper branch, the 
high n+, the specific heat is negative, and as the 
heating rate and the luminosity increases, the 



temperature decreases; in the lower branch, the low n+, 
as the heating rate increases, the temperature increases 
and eventually reaches a T c , where both branches 
merge. For a further increase in the heating rate, the 
behavior follows that of the high n + branch, with 
temperature decreasing till the plasma becomes 
effectively thick. These findings of Lightman[8] have 
been confirmed independently by Svensson 1982[17], 
1984[18]. 

When these results are applied to flows in accretion 
disks, the role of the temperature on the existence of 
equilibrium production-annihilation is now played by 
the accretion rate. If flows in accretion disks have 
conditions such that the dynamical time is shorter than 
the characteristic time for the ions to transfer energy to 
the electrons, radiative cooling will be inefficient and 
most of the energy will be advected with the ions. A 
two-temperature disk will follow, with both T ; and T e 
close to their virial values, with T ; » T e ~m e c 2 , close 
to the inner radius. One, then, should expect a 
significant pair production in these flows. Kusunose 
and Takahara 1988 [9] have shown that, for accretion 
rates below a critical one, there exist two branches of 
solution under equilibrium production-annihilation: 
one of them, the high pair, has z>l, and the other, the 
low pair solution, has z<l. For accretion rates greater 
than the critical one, no equilibrium solution is 
possible. An interesting result from these works is that 
the existence of a critical accretion rate is independent 
of pairs being created by particle-particle, particle- 
photon, or photon-photon processes, as long as the 
photons are internally produced. When photons are 
externally produced, the rate of pair creation is no 
longer non linear in the particle density. 
As far as important issues like the number of solutions 
accessible to the disk for given physical conditions, 
the role of the accretion rate and viscosity in 
determining steady state solutions for accretion disks 
with pairs produced by internal photons are concerned, 
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reasonable agreement is achieved by several 
authors(citation). However, if one is concerned with 
apir production by external photons, the situation is 
less clear, and the results are even discrepant. Tritz 
and Tsuruta 1989[19] and White and Lightman 
1989[21] obtain rather similar results for disks with 
pairs produced by external soft photons. According to 
them, pairs produced in that way do not substantially 
modify the disk, the density of pairs being negligibly 
small. Since they failed to find a critical accretion rate, 
pair production-annihilation equilibrium is always 
possible. Their solution are not self-consistent in the 
sense that the ion temperature exceeds its virial value 
and the disk scale height is larger than the radial 
distance. Kusunose and Takahara 1990[10], on the 
other hand, relaxing the imposition of setting the y 
Kompaneetz parameter equal to one, have studied 
disks with pair production by external soft photons, 
with different soft photons input. They found a crtical 
accretion rate and pair densities smaller than the 
proton density. The critical accretion rate is 
proportional to the viscosity parameter and depends on 
the soft photon energy. 

They also have found a unique solution to the disk 
equations. However, their result is only numerical and 
they have exploited only a small range for the soft 
photon energy, and they have not included advection 
in their energy equation. The effect of advection was 
later included by Kusunose and Mineshige 1996[11]. 
They, however, considered only internal soft photons 
produced by bremsstrahlung. Misra nad Melia 
1995 [12] have also considered pair production by soft 
external soft photons and internal bremsstrahlung 
photons. They have included a term of energy cooling 
due to photon escape, but they have not included 
advection. 

In recent years, a large amount of observational data 
on compact X-ray and y-ray sources are highly 
indicative of electron temperature close to the electron 
rest mass and ion temperature close to its virial value 
(citation) and, in the Inter Stellar Medium, close to 
some sources, there evidences of pair annihilation 
features (citation). This, besides being suggestive of 
the presence of pairs, implies not the importance of 
advective cooling, due to the increase of the scale 
height as compared to the radial distance, but as well 
that the inner region, where the disk has thickened, 
will be subject to a large amount of soft photons, 
irradiated from the outer disk. 

If the picture we have about these sources, surrounded 
by accretion disks, is correct, the irradiation of the 
inner region is unavoidable. Therefore, it will be very 
useful to revisit disks with pairs produced by external 
soft photons. We, then, propose, in this paper, to 
undertake this task. Our special emphasis will be: 

The search for multiple solutions for soft 
photon comptonized two temperature 



accretion disks, with photons supplied by an 
external source; 

The search for critical values for both the 
accretion rate and the viscosity parameter; 
The search for conditions under which the 
solutions are physically consistent; 
The search for maximum electron 
temperature at each radius. 
The outline of this paper is as follows. In§ I we 
present our model of accretion disk, together with 
some assumptions. In §11 we explain the units and 
obtain the disk equations. In §111 we obtain the 
solution to the disk equations. In § IV we obtain the 
critical parameters. In § V the conclusions are 
presented. 

LABOUT THE MODEL 

The model of accretion disk we shall be considering is 
a soft photon Comptonized two-temperature accretion 
disk, with pair production by photon-photon 
interaction. The photons are externally produced and, 
after impinging the disk, they are upscattered in 
energy. The region we are most interested is the 
innermost one, where the disk has thickened. Most of 
the standard assumptions will be kept: keplerian 
velocity, stress tensor ~aP, a being the viscosity 
parameter and P the pressure, hydrostatic equilibrium 
in z-direction, particles obey a Maxwell-Boltzmann 
distribution. The energy transport will be radiative in 
z-direction, perpendicular to the plane of the disk, and 
advective in the radial direction. The advection of 
energy will be treated locally. The coupling between 
protons and electrons will only through Coulomb 
collisional energy exchange. 

For the radiation field, we shall follow the treatment of 
Svensson 1984[18], and Zdziarski 1985[20], and 
assume it may be represented by an inverse power law 
with an exponential cutoff at KT e plus a Wien bump. 
For the lepton energy equation we shall follow the 
Comptonization treatment as given by Zdziarski 
1985 [20], which is equivalent to to obtain a relation 
between the spectral index y, the electron temperature 
and the probability of a power law photon into the 
region x(=hv/ KT e )>9 e (KT e /n^c 2 ), where the Wien 
spectrum is formed. 

As usual, it is assumed that heating of the disk by 
viscous processes locally balances cooling. 

II.DISK EQUATIONS 

Throughout this article we shall adopt the following 
set of units: radial distance and the disk scale height 
are expressed in units of R s =2GM/c 2 , where R s is the 
Schwarzschild radius; electron temperature 8 e , and 
proton temperature 9 P in units of m e c 2 ; M is the mass 
of the central object in units of 10 solar masses;ltl is 



the accretion rate in units of L E /c , where L E is the 
Eddington luminosity; z, the pair density will be 
expressed,in units of the electron number density N e . 
In the following we shall adopt the same procedure of 
Bjornsson and Svensson 1992[3], and Bjornsson et al. 
1996[4], and employ, as much as possible, as 
variables, the pressure P in units of rest mass energy 
density, the compactness parameter I , or radiative 
flux in units of 47im e c 3 /(3cTH), and the Thomson 
scattering depth, x p = o T N e H, where a T is the Thomson 
cross section for electron scattering, and H is the disk 
scale height. 

Using these variables, the disk structure equations, 
hydrostatic equilibrium, conservation of angular 
momentum and energy balance are written, 
respectively, as 
1H 2 

P = - — , (1) 



2r 3 ' 



(2) 



Now, defining r\ in the same way as Bjornsson and 
Svensson 1996, i.e., 

rj = □r" 15 5(r), (4) 
And substituting the scale height by the pressure, we 
obtain 

i v 



And 

□ - -. — i4nP , 

V2m ' 

(6) 

Where 

^ Qadv 

Q+ ' 



(5) 
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Q + , and Q a(iv being, respectively, the flux of heat 
generated by viscous processes, and the advective 
flux. 

For the advective cooling, we shall use an approximate 
local expression (Abramowicz et al. 1996[1]), 

" 2 (8) 



Qadv = Q+(.-r- 



Therefore, 
A — 1 — 2Pr. 



(9) 



It should be remarked that eq.(7) is equivalent to 

F r = Q + -Qadv, (10) 

Where F r is the radiative cooling. In steady state, the 
leptons emit exactly the energy they receive from the 
protons through Coulomb collisional energy exchange, 
which may be written, approximately, as (Stepney and 
Guilbert 1981 [14]) 



F r = 4.54 x 10 26 j^p (0 p - e )(l + 2z) ( 1 + 
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The energy equation for electrons and pairs comes 
from the treatment of Comptonization, as given 
Zdziarski 1985, and it reads 

1 + 8° 5 

T = a + 2z)T * = Wa^W) ■ (12) 

In writing eq.(12), we have assumed y=l, which is not 
very far from observed spectral indexes in a large 
number of X-ray and y-ray sources. Substituting 
eq.(12) into eq.(ll), we obtain from the definition of 
the compactness parameter 



4n j H 
t = — F 

3 m e c s 



(13) 



(i + o^y 



* = °oso4T p (e p -e e ) -; e m }y (14) 

Now, writing the equation of state for a gas composed 



(3) of protons, electrons, pairs and radiation, i.e. 



m„ 



p =— (e p + (i + 2z)e e +— (i + t w ), 



TO. 



4ti 



(15) 



Where m e ,m p are, respectively, electron and proton 
masses. x w is the Wien averaged scattering optical 
depth (Svensson 1984[18]), and it reads 



(1 + 0° 5 ) (1 + S6 e + O.40 e 2 ) 



(1 + 1202) 
(ln(1.126» e ) + 0.75) 



e e <i 



e R > i. 



(16) 



U6 e 25 (l + e )(l + 120 2 )' 
Before we write down the equation for pair 
production-annihilation equilibrium, we must talk 
photon balance, from which we can have the radiation 
field to obtain the rate of pair production by photon- 
photon interations. 

Following Svensson 1984[18], we treat the radiation 
field as 

N = e-*(^x- 2 +^x 2 ), (17) 

Where we have specialized for an spectral index of 1. 
The ratio of intensities is given by Sunyaev and 
Titarchuk 1980[16], i.e., 

~Np~ ~ r(5) Psct ' (18) 

Where F is the Gamma function, P sct is the probability 
of scattering, and N p is related to the soft input N soft by 
conservation of the total number of photons. 
If N SO ft is the rate of soft photon creation, and we have 
equilibrium between input of soft photons and output 
ot photons from the disk, 

□ y = n so/t , (19) 
And 
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H 



(20) 



Since the photons escape from the disk with an 
energy~m e c 2 e , 
Fr 

N v = —(1+t). (21) 



7 mec 3 0e 
However, 

N y = I N(x) dx, 
Jx 



(22) 



Where x is the soft photon energy in units of kT e 
.Therefore, 



x A e x dx + ^Y{^) T 



F r 



(1+T), 



(23) 



N„ = 



mec 3 9e 

Where we have used (Zdziarski 1985[20]), 

Psct = T. (24) 
Now, using eqs.(13), (14), (21), and (24), we may 
write 

x~ 2 e~ x dx + 0.258^ 

x-(1+t), (25) 

Which gives for the pair creation-annihilation 
equilibrium equation (Lightman 1982[8]; Svensson 
1982[17]) 

(1+t) 2 ^/(x o ,0) = B(9)(t 2 -t 2 ), (26) 
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An 2 
Where 

B{G) = ■ 



In (1.126 + 1.3) 



(27) 



20 2 +n (1.126+1.3) ' 

f(x ,6) = 

0.693 (S™x~ z e-*dx0.258T) 

[8 S °+ (0.128t) 2 Q~ x + 0.128t(6» + 6> 3 )]. (28) 
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III.THE NUMBER OF CONSISTENT 
SOLUTIONS TO THE DISK EQUATION 

Using eqs.(5), (6), (14) and (15), we finally reduce 
our system of equations to 

Q5 = 4~(l+g e 05 ) 2 
V2m e V (1 + 1201) 

(1819.2Pt p - T p 6>)(l + 120|) - 9° 5 (1 + 9° s ) 
(1 + t w )(1 + 6» e - 5 ) 2 + 249.336> 2 (1 + 126> e 2 ) 



(29) 



And 
9 m 



Where A, r\, x, f(x ,9), B(9), x p are given, respectively, 

byeqs.(4),(7),(12),(15),(26),(27). 

For given T|, a, x ,r, this is a system of two-coupled 

equations on the variables P and 9. 

We now solve this system of equations specializing 

for Cygnus X-l, for which we assume canonical 

values of M=10M and m=6.52xlO" 2 (Liang and Nolan 

1984[7]). We, then, obtain the following solutions for 

the different sets of parameters. 







P 




0.621 


34.34 


0.0195 


0.239 


1.201 


171.56 


0.0958 


0.0485 


z 


A 


H/r 


i 


0.345 


0.805 


0.441 


0.597 


0.576 


0.042 


0.979 


0.0693 



Table(l). The solutions to the disk equations for the 
parameters r=5, rr=1.315x 10" 3 ,oH).l, x =0.1 







P 


Xp 


0.715 


29.185 


0.0121 


0.181 


1.085 


353.07 


0.0300 


0.073 j 


z 


A 


H/r 


I 


0.345 


0.879 


0.348 


0.243 


0.385 


0.700 


0.547 


0.304 



Table(2). The solutions to the disk equations for the 
parameters r=15, r|=6.206x 10" 4 ,a=0.1, x =0.1 



©e 


©p 


P 


Xp 


0.621 


34.166 


0.0194 


0.240 


1.365 


84.46 


0.0987 


0.0471 


z 


A 


H/r 


I 


0.342 


0.806 


0.440 


0.597 


0.344 


0.013 


0.994 


0.0231 



Table(3). The solutions to the disk equations for the 
parameters r=5, rr=1.315x 10" 3 ,aO.l, x =0.01 



©e 


©p 


P 


x P 


0.715 


29.03 


0.01207 


0.182 


1.13 


420.36 


0.03107 


0.0706 


z 


A 


H/r 


I 


0.341 


0.879 


0.347 


0.242 


0.342 


0.689 


0.557 


0.305 



Table(4). The solutions to the disk equations for the 
parameters r=15, r)=6.206x 10" 4 ,oH).l, x =0.01 



©e 


©p 


P 


Xp 


0.546 


13.10 


0.00761 


0.305 


1.85 


1518.3 


0.0956 


0.0243 


z 


A 


H/r 


I 


0.341 


0.924 


0.276 


0.428 


0.348 


0.044 


0.978 


0.0719 



Table(5). The solutions to the disk equations for the 
parameters r=5, r)=1.315x 10" 3 ,oH).2, x =0.01 
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®„ 


P 




0.597 


7.79 


0.00425 


0.258 


1.603 


2894 


0.0329 


0.0334 


z 


A 


H/r 


i 


0.341 


0.958 


0.206 


0.156 


0.343 


0.671 


0.574 


0.305 



Table(6). The solutions to the disk equations for the 
parameters r= 15, r)=6.206x 10" 4 ,a=0.2, x„=0.01 

The first thing to be remarked is the existence of two 
solutions. For both them the number of pairs does not 
outnumber the protons. These solutions are, 
practically, independent of the energy of the soft 
photons. Depending on the combination of the 
parameters the solutions are physically consistent 
close to the hole. For a=0.1, every solution we found 
is consistent for r=5. If we go farther out, after a 
certain radius, the branch which has the greater 
temperature will become inconsistent, i.e., the proton 
temperature will be greater than its virial value, or the 
disk scale height will be greater than the radial 
distance. We haven't found this radius, but certainly it 
is inside r=15. For a=0.2 there is only one physically 
consistent solution: the branch with the greater 
temperature does not meet consistency. The variables 
have quite different behaviors in different branches. In 
the lower branch (the one with the smaller 
temperature), the electronic temperature decreases and 
the proton temperature increases as we come close to 
the hole. In the lower branch, the pressure, the scale 
height, the Thomson scattering depth, and the 
compactness increase as we go into the inner region; 
the radiative efficiency decreases, and the pair number 
is practically constant. In the upper branch, the 
pressure, the scale height, the pair number increase as 
we go closer to the inner radius; the radiative 
efficiency, the Thomson scattering depth, and the 
compactness decrease. 



IV.CRITICAL CONDITIONS 



We now look for critical conditions for the disk, i.e., 
parameters values which define the transition from a 
steady state to a time dependent solution. In order to 
do so, we start with eq.(29) and treat it as a function of 
P, assuming all other variables as parameters. 
Besides, a cursory analysis reveals that the term 
x p 9 e (l + 120|) is negligible in most situations. We, 
then, write eq.(29) as 
3.22 x 10 2 r)AP 0S = g 0l 
Where 



(31) 



do = 



643.18^(1 + 120|) + 0° 5 (1 + 0° 5 ) 
(1 + t w )(1 + 0°- 5 ) 2 + 249.330|(l + 120|J 



(32) 

The left hand side of eq.(31) has a maximum at 
P=l/6r, which equals 87.8/ Vr . Then, in order to have 
solution , this maximum should be greater or equal to 
the right hand side of eq.(31), i.e., 

> go- (33) 
Then, criticality is determined by the value of a: if a is 
greater than a c defined by the equation below 

a = 7.35Vr(l+e°- 5 ) 2 

C (l+T W )(l + e^ 5 ) 2 +249.338|(l + 12e|)' 

If a is greater than a c , the system of equations for the 
disk will admit two steady solutions, no matter the 
value of T|. However, if a<a c , the system has a critical 
condition, and it admits two solutions only and if only 
r|i< T|<r| 2 , with 

0°- 5 (l + 0°- 5 )Wr 
Vi = — r- 2 < (35) 



9i 



And 

V 2 = 
where 



2.0905(0)° Wr 

32 



(36) 



g 1 = 7.35(1 + 0° 5 ) 2 (1 + 120|)Vr - a (1 + 120|) 
X [(1 + T W )(1 + e 05 ) 2 + 249.330 e 2 (l + 120 e 2 )], 

eq.(37) 

And 



32 = 



19.739— + 

a 2 



I 2 1.21 xlO 6 



(1 + t) 2 /(x o ,0) 



(38) 
let us 



To better understand these critical conditions, 
take r=5. 

For this value of r, we see that the maximum value a c 
can attain is ~2. Therefore, for any a smaller than 2, 
solutions will be steady outside a region in 9, which 
borders are determined by the value of a. Let us 
specialize for <x=0.1. The critical region is 3. 
7954xl0" 5 <9<0.612. We now ask if there is a 
solution for r|=0.01 and 9=0.5. For these values of 
a,9, and r, we have r| outside the range r|i=1.02x 10~ 4 
and ru=5.69x 10" 3 . Using eqs.(4), (35), and (36) we 
may write 



< n < 



0»- 5 (l + e a5 ) 3 ar 2 



5l S(0) 
V.CONCLUSIONS 



2.O90B(0) O5 ar 2 



g 2 S{B) 



(39) 



The results we have obtained differ from those of 
White and Lightman 1989[21], Kusunose and 
Takahara 1990[10], and Kusunose and Mineshige 
1996[1 1], mainly by the number of solutions to the 
disk equations and by the conditions under which the 



system reaches criticality. The critical viscosity 
parameter a c , below which the system may develop 
time behavior, is a function of r and 9 e . For a given a, 
we may find a region in 9 e where a<a c . Even under 
this condition, the system may have two steady state 
solutions, as long as the accretion rate satisfies ltl| <ltl 
<m 2 , where ltl| and ltl 2 are functions of a, r, and 9 e . To 
make these points more clear, let us take a=0.1 and 
i=5. We, then, obtain <x c =0.1 at 9 e =3.8xl0~ 5 and 8 e 
=0.612. For 3.8xl0" 5 <8 e <0.612, a c >0.1. Therefore, 
for 8 e =0.5, and a=0.1, there will be two steady 
solutions only and if only 5.06xl0" 3 <rtl <0.282. If we 
take, at r=5, values of 8 e , in that range, that maximize 
a c ,rfli, and minimize ltl 2 , we obtain a c ~2.0, 
rhi~0.033, and ltl 2 ~0.2. As compared to White and 
Lightman 1989 [21] results, ours reflect the fact we 
have included advection in the energy equation and 
have adopted a different approach to the treatment of 
Comptonization. It seems that the inclusion of 
advection in the pair equation, as White and Lightman 
1989[21] did, is a minor effect as compared to the 
inclusion of this term in the energy equation. We have 
included advection in a rather local manner 
(Abramowicz et al. 1996[1]), explicitly dependent on 
r, and this has introduced an explicit r dependence in 
our equations, and not only an implicit dependence 
through n. This somehow helps to explain 
Our different result as compared to Kusunose and 
Mineshige 1996[11] results, besides the fact they have 
used different radiative cooling, and used another local 
approach to the advection treatment. Though 
Kusunose and Takahara 1990[10] have not included 
advection in their disk with pairs produced by soft 
external comptonized photons, they obtain a rather 
similar result to the maximum accretion rate. Their 
results are, however, very sensitive to the soft photon 
energy. Our results qualitatively ,are in agreement 
with those of Chen et al. 1995[5] ,who have analyzed 
solutions to the disk equations in terms of itl and a 
and have found a critical a c , above which there is 
always at least one optically thin advection dominated 
solution for a given accretion rate. For a<a c , there is 
a maximum accretion rate, above which there are no 
steady state solutions. We, finally, may summarize the 
main results obtained in this work, in which we have 
revisited the inner region of a two temperature soft 
photon comptonized accretion disk with pairs 
produced by photon-photon interactions: 

We have found two solutions to the disk 

equations; 

For both these solutions the number of pairs 
is always smaller than the number of protons; 
These solutions practically do not depend on 
the soft photon energy; 
The consistency of the solutions close to the 
hole depend on the parameters combinations; 



Going farther out, after a certain radius, the 
upper solution is no longer consistent, i.e., the 
proton temperature will be greater than its 
virial value, or the disk scale height will be 
larger than the radial distance; 
In the lower branch, the one with the smaller 
temperature, the electronic temperature, the 
radiative efficiency decrease, the pair number 
is practically constant, and the proton 
temperature, the pressure, the scale height, 
the scattering depth, and the com pactness 
increase as we come closer to the hole; 
In the upper branch, the electronic 
temperature, the pressure, the scale height, 
the pair number increase ; the proton 
temperature, the radiative efficiency, the 
Thomson scattering depth, and the 
compactness decrease as we come closer to 
the inner radius; 

If a increases, temperature decreases in the 

lower branch and increases in the upper one; 

If m increases, temperature remains 

practically constant in the lower branch and 

increases in the upper one 

For a=0.2, only the lower branch solution is 

consistent; 

In the lower branch, cooling is dominated by 
radiation and, in the upper branch, cooling is 
dominated by advection; 
It may happen that at certain r, both solutions 
are consistent. When we go farther out, just 
one of these solutions keeps consistency. 
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